Prepared for submission to JCAP 



Inflation and primordial non-Gaussianities 
o|of ''generalized Galileons" 

o 

CD 

Q 

g" Xian Gao"^^ Daniele A. Steer 

Q_J \ "Astroparticule et Cosmologie (APC), 

_7 ■ UMR 7164-CNRS, Universite Denis Diderot-Paris 7, 10 rue Alice Domon et Leonie Duquet, 75205 Paris, France 

O . ''Laboratoire de Physique Theorique, Ecole Normale Superieure (LPTENS), 

' ' 24 rue Lhomond, 75231 Paris, France 

^ . '^Institut d' Astrophysique de Paris (lAP), 

'^ ' UMR 7095-CNRS, Universite PieiTe et Marie Curie-Paris 6, 98bis Boulevard Arago, 75014 Paris, France. 

cd :' 

I— —I, E-mail: xgao@apc.univ-paris7.fr, steer@apc.univ-paris7.fr 

►^ , Abstract. We set up cosmological perturbation theory and study the cosmological implications of the so-called "generalized 

pg ' Galileon" developed in [6, 7]. This is the most general scalar field theory whose Lagrangian contains derivatives up to second 

T^ , order while keeping second order equations of motion, and contains as sub-cases fc-inflation, G-inflation and many other models. 

^^ We calculate the power spectrum of the primordial curvature perturbation, finding a modification of the usual consistency relation 

^^ . of the tensor-to-scalar ratio in fc-inflation or perfect fluid models. Finally we also calculate the bispectrum, which contains no 

C — ' new shapes beyond those of fc-inflation. 

o 



% 



Keyword SI modified gravity, inflation, physics of the early universe 

ArXiv ePrint: 1 107.2642 



Contents 

1 Introduction 1 

2 Generalized Galileons: Field equations 2 

3 Background equations 4 

4 Linear perturbations 6 

4. 1 Power spectra 6 

4.1.1 Scalar perturbations 6 

4. 1 .2 Tensor perturbations 7 

4.1.3 Quantization 8 

5 Bispectrum of curvature perturbation ( 9 

5.1 Cubic action for the curvature perturbation 9 

5.2 Bispectrum of the curvature perturbation 12 

6 Conclusion 13 
A Equations of motion 14 
B Explicit expression for T^i/^ 15 
C General consideration of the structure of perturbative action 15 

1 Introduction 

Recently a great deal of effort has been devoted to developing consistent modified gravity theories which, as an alternative to 
to dark energy or the cosmological constant, may provide an explanation of cosmic acceleration (for a recent review, see [1]). 
Such IR modification of gravity arise, for example, in higher-dimensional setups as in the DGP model where self-acceleration 
is sourced by a scalar field (j), the helicity-0 mode of the 5D graviton. On small scales the DGP model also reproduces general 
relativity due to non-linear interactions operating through the Vainstein mechanism [2]. Other theories of modified gravity, such 
as DBl-galileons can also be obtained from a higher dimensional approach [3]. 

Generally, modified gravity models contain additional degrees of freedom which can often be viewed as scalar fields. In 
this paper we consider models with a single scalar field (f>, though allowing for couplings between (j) and the metric which may 
be very different from the standard couplings and potentials considered in many inflationary models. Indeed, the context of our 
work is that of the "generalized Galileon" derived in [6] (see also [7]). 

While there are numerous different models of inflation with different actions and potentials [8] - for example chaotic 
inflation, small/large-field inflation, /j-inflation and DBl-inflation to name a few - single-field fc-inflation [9], in particular, at- 
tracted a lot of attention since it is the most general scalar field theory with a Lagrangian containing derivatives up to first order 
£ = €{(/), V(/>) (and hence, clearly, having second order equations of motion). Indeed, all the inflationary models listed above 
are a special case of fc-inflation, and for that reason the development of cosmological perturbation theory in general case of 
fc-inflation has a broad range of applicability. 

In this paper we go beyond fc-inflation and study the most general Lagrangian in 4 space-time dimensions which contain 
both first and second derivatives of the scalar field, C = C{(j), V(/), VV(/)), but constructed such a way that the field equations for 
both the metric and the scalar field remain second order [6, 7]. This prevents the theory from having extra degrees of freedom as 
well as Ostrogradski instabilities [10], and thus yields a possibly viable higher-order derivative scalar field theory. This general 
Lagrangian includes, in certain limits, the decoupling limit of DGP model [11] as well as some consistent theories of massive 
gravity [12, 13]; and it also includes the "Galileon" model [14], the conformal Galileon [15], K-mouflage [16], and also "G-" or 
"KGB" inflation [20, 33]. Finally, in a very simple limit (which can be useful to check results) it also reduces to fc-inflation. Our 
aim is to study cosmological perturbation theory in the context of this very general Lagrangian, outlined in Section 2, working 
both to second as well as third order so as to study non-Gaussianities. The formalism we develop is therefore applicable to a 
large class of models. As such, this paper is rather technical - since we crank the standard handle of cosmological perturbation 
theory - but the results are rather general. 



As we will see below, see equation (2.1), the Lagrangian we consider [6] contains four free functions which will play an 
important role in the following. We denote them by 

K{X,(t)), and G("'(X,0), n = l,2,3 (1.1) 

where 6 is the scalar field, and 



X^--g^"'d^^d.cj), (1.2) 



1 
2 
is its kinetic term (we work with a mostly positive signature). To obtain the Lagrangian of fc-inflation [9] one must simply set 



K^O, G(^'2'^)(X,0) =0, (fc-inflation limit). (1.3) 

The Galileon Lagrangian [14, 17] is obtained when 

K = Xc(o), G'-^^ = Xc(i), G^^/^ = Xc(2,3), (Galileon limit) (1.4) 

where the C(„) are (dimensionful) constants. (The effective scalar field Lagrangian in the decoupling limit of the DGP model is 
obtained by setting C(2) = C(3) = 0.) In this case the Lagrangian will be invariant under shifts of the field, namely </>—)- + c. 
In general, when K{X, (j)) and G^^'^X, (f>) are non-vanishing functions of X and 0, shift symmetry is broken (as indeed is the 
"Galileon" symmetry [14] since we work in curved backgrounds). Without shift symmetry, questions of the stability of the theory 
against large renormalisation may arise [28], but we do not dwell on this point here. 

To obtain the Lagrangian of G-inflation [20] (or equivalently KGB-inflation [33]), one must set 

K,G^^^ y^ 0, G(2) = G^^) = (G-inflation limit). (1.5) 

Cosmological perturbations has been studied in Galileon model [ 1 8-25] . Primordial non-Gaussianities [26-3 1 ] have been studied 
in G-inflation limit, and it was shown [29] that to leading order in slow-roll parameters there are no new momentum shapes for 
the bispectrum (beyond the standard ones of fc-inflation). Despite that, a distinctive feature of G-inflation is that, though it has 
only a single scalar field, its energy momentum tensor takes the form of an imperfect fluid [33, 34] (contrary to the fc-essence 
models). This fact sheds some light on the understanding of imperfect fluids in cosmology, and furthermore can lead to violations 
of the null enery condition thus opening up a whole range of possibly exotic phenomena [35] . Problems of superluminality may 
arise in this case, though it is still under debate as to whether causality is violated in general [36]. Previous understanding of the 
conservation of curvature perturbation on super-Hubble scales relied highly on fc-essence or perfect fluid models: it is interesting 
to see [38, 39] that this conservation law still holds in the generalized Galileon model, that is when the functions K and G^"^ are 
all taken to be arbitrary functions of (j) and X. 

In this work, we keep the greatest generality by taking the full Lagrangian of generalized Galileon [6, 7] into account: that 
is, arbitrary functions K and G*^"-* . We derive the background equations of motion for an inflationary background, around which 
we calculate the power spectrum and the bispectrum of the primordial curvature perturbation C, (the background evolution and 
Unear perturbations of the generalized Galileon [6, 7] were also investigated very recently in [23], see also [24] for a discussion 
of self-tuning mechanism on FRW background). We find a modification of the familiar consistency relation of the tensor-to- 
scalar ratio. From the theoretical point of view, the results presented here can be used to study the evolution of cosmological 
perturbations as well as the observational impUcations numerous models. 

This paper is organized as follows. In Sec. 2, we review the generalized Galileon model. Background evolution of our model 
and linear perturbations are studied in Sec. 4. In Sec. 5, we compute the full third-order action for the curvature perturbation and 
evaluate the corresponding contributions to the bispectrum. Final section is devoted to conclusion and discussion. 

2 Generalized Galileons: Field equations 

In [6], the most general scalar field action which contains at most second derivatives of (p, is polynomial in these second deriva- 
tives, and which leads to second order equations of motion in flat space-time, was constructed in D dimensional space-time. Its 
"covariatization" [17] leads to an action with field equations which are second order or lower for both and the metric. In 4D 
this is the unique action with second order equations of motion [7], and is given by 

(2.1) 




where we have set Mp, = l/(87rG') = 1, and 
£(1) = g(^) (X, 0) D,/), 



;n0)2-(v^v,0)' 



£(3) ^ g|^'(x» (n0)3-3n0(V^V,,; 






(2.2) 
(2.3) 
(2.4) 

(2.5) 



Here, as discussed in the introduction, K{X, cp) and G^"^ (X, 0) are arbitrary functions of (j) and X (defined in (1.2)). In C'-^K 
(V^Viy0) = [^ fj^ u4>) (V^'Vc/)) and i? is the Ricci scalar. Thus in principle the Einstein-Hilbert action R/2 in (2.1) can be 
viewed as a special case of £(^^ with G*-^' = 1/2. However, since for certain applications of our analysis we may wish to set 
£(2) _ Q (^ggg (1.3)-(1.5)), we keep the Einstein-Hilbert term explicit. Whilst the Lagrangian £'°) and C^^^ axe. familiar from k 
and G-inflation, terms C^'^'^^ represent completely new class of scalar field theories which has not been explored so far In C^^\ 



(V,,V^<i 



{V^N^,cp) (V^V^) (VpV^ 



whereas G^u is the Einstein tensor. Notice that a term Gf^^d'^cfid^cl), discussed in for example [25], can by obtained by setting 
G(3) ex (/). 

For completeness, we collect the corresponding equations of motion for the Galileon scalar field and the metric. The 
Einstein equation following from (2.1) is 



^ nu ^11 






n=0 



where a superscript "(")" denotes the contribution from £^"). Using the notation G^ = dG'^"'^ /dX and G 
stress-energy tensors TJ),^' ' are given by [33, 34] 



.(«) 






(2.6) 
aG("V5<A.the 

(2.7) 



G^yVa^V^ 



2XGy ) 5m, + (G^,ljn^ + 2G[IA V^0V,0 + 2GWV(^0V,)X. 



^(1), 



(2.8) 



Here, from (1.2), V^X ~ — V^Va^V 0. As observed in [33, 34], Tf^J contains second derivatives (but not higher, by construc- 
tion) of (j), and as a result Einstein's equations (2.6) (as well as the scalar field equations of motion) contain second derivatives 
of both the metric g^„ and of (f) as soon as £(^^ is present. There is no conformal transformation that can diagonalise the system 
of equations, which remain coupled: this phenomenon has been dubbed kinetic braiding [33]. The expressions for Tp^ii are 
rather more involved, again containing second derivatives of (p, but they now also contain second derivatives of g^j.^. In particular 
[23, 39] 

TlS = <?M-{i?G(2) - G J) ((00)2 - (V,V,0)2) + 4XG(^^ - 2G%ypX^PX 

^2 i2G% + D^Ggx) V.^VX + 200 (2XG% - gJJ) + 2G;^'i?p,V''0V-0} 



G^^^i? + 4n0G*2^ 



,X0 



2G 



(2) 



G 



(2) 
XX 



((00)2 - (VpV„ 



Vu0Vt.(5 



+4 (o0gJ^^_:^ + 2gJxV) V(a'</'V-)^ + 2G^i'x (V^XV.X - 2VpXV''V(^0V,) 

+2 (g,'^'O0 + G|l.';^Vp0V''X - 2XgJ2^^ + gJ ) VpV,0 

-2G|^^ (V''Vp0V,Vp0 + 2V(p0i?,)pV''0 + i?p^<,,V''0V"0) - 2G(2)ii-^,. 



(2.9) 



.(3) 



The expression for T^^J is given in Appendix B. 

While the stress-energy tensor T^,^ is clearly of the perfect fluid form, T^,' ' ' are not. In the local rest frame defined by 
the effective four- velocity 



/2X 



Uf.U^' = -1 



(2.10) 



the stress energy tensors take the general form 

T^:^ = P^'^^u^u. + P^'^Hg^. + u,u,) + (w^gl") + «.9,^r)) + ^("j (2.11) 

where p^"'' = Tj^i, u^u^ is the energy density, P*^") the isotropic pressure, qjl^ the energy flow, and ttviI = tt/j" the anisotropic 
stress. The expressions for p'^^\ P^^'^ and (/'^^ can straightforwardly be read off from (2.8) from which it follows that tt^,; = 0. 
For Tliv the situation is yet more complicated, and from (2.9) and (B.l) they both lead to non-vanishing anisotropic stress 
■Kiiv as well as energy flow. The consequences of this non-perfect fluid picture have been discussed in [33, 34, 39]. 
The equation of motion for the scalar field (f) is given by 

3 



n=0 



0, 



where 



£("' =v„j("'^ + /: 



in) 



The currents J^^\ which are only conserved when the theory is shift-symmetric CV = 0, are given by 



j(o)m 



,x' 



j(i)^ = n(/)G4iv^0 + V^G(i), 



j(2)m ^ 

j(3)m ^ _ 



.(2) 



.(2) 



[Uc^Y - (V,,V,0)^ G\i>^ + RG 



G<^) (3 [Ucjyf g^- - -ig^- {Vp^Af - 6^0V^V^ 



2V^(G^^^(n0g^''-V^V''(j 



G 



(3) 



Av7l', 



GV^Va^V-'V 



-GG^'^V.G^^^ 



(2.12) 
(2.13) 

(2.14) 
(2.15) 
(2.16) 



(2.17) 



Though at first sight it might appear that these equations of motion contain derivatives of order 3 or higher, this is not the case: by 
construction, these higher order derivatives cancel on calculating V^ J^"-'^. The explicit equations of motion given in Appendix 
A. 



3 Background equations 

We now consider a spatially flat FLRW geometry with metric 

ds2=:a2(r/)(-dr/2 + da;2^, 

where 77 is conformal time. Then the kinetic term X reduces to 

1 



--^^2^*° 



(3.1) 



(3.2) 



where a dash denotes a derivative with respect to r]. Then the background equations of motion (which most straightfowardly are 
obtained by the requiring first order perturbative action to be vanishing to vanish) are given by 



and 



^ =y^ 



-H^ + 2-H' + a'P = 0, 



(3.3) 
(3.4) 



where the effective energy density p and pressure P can be read off from (2.7-2.9) as well as (B.l) and are given by 



n=0 



P = ^P("), 



n=0 



(3.5) 



-4- 



with 



3(°' = 2XK ^-K 



X 



p(i) = _6^XGW+2XG 



(1) 



■,(2) 



.n' 



6lL G^2) _ 4^ ^(2) ^ xG^'J^] - 6^ (2XG 



,x 



.w 



.(2) 



,'H 



(3)=36^X 2XG(1 + 3G(!) 



a^ V 



,X0 



12HV 
a4 



G 



(2) 



,X</) ^ ^,ct> 



X ( 5G|^) + 2XG(^);, , , 



and 



P(i) = 2XGy + ^^gJ (0" - H0') , 

P(2) = 4XG Ji + ^G(2) (^2 + 2H') + ^XG% (0" - 3H<^') + ^G^^ {H^' + </)") 



tHXG^^J^ [2HX 






p(3) 



rJiXG 



,xxy---- „2 
24-^2 



(3) A' 



24, 

'~2' 

-^HXG|i^^ ( 3HX 



2^,X 



AX (n^ - 2H') - ^(jy'q 



X'^G)-^xx ((('" - '^(f'') + -^G'^x^ [(2'H' - S-H^) 0' + 3H(t)"] 






(3) 



X (3H'' - 2H') - ^</)>" 



(3.6) 
(3.7) 

(3.8) 

(3.9) 

(3.10) 
(3.11) 

(3.12) 
(3.13) 



Of course the energy fluxes and anisotropic stresses vanish on the background, and (3.6)-(3. 13) reduce to the standard expressions 
[9, 14, 17, 20] in the fc-essence (1.3), Galileon (1.4) and G-inflation limits (1.5) respectively. 

The equation of motion for the scalar field follows from (2. 12) (or alternatively as a combination of (3.3) and (3.4)) and is 
given by 

— (aV)' = £,0 



with 



J = K,f-JpGW+3HG(iV)-^ 



•HG' 



(2) , A' fM2) r<(2) '^<^' 



X 



■-^ X<p '-'.XX ^2 



en 



2aI r 



-ia^G'-l^ + (j)' { 3a^'HG[^J + cj)' (-Sa^G*'^^^ + nG[^J^(j}' 



and 



C 



4 = a;^ - IgJ {2h<p' + 0") + le (a^d^ {n' + W) + nG%<i>'r) 



n6 



HG^t^^ ((/)')' {-2W + 30") + Sa^G^'^ (2H>' + n4>") 



,Xd 



As usual, we define the slow-roll parameter 



d\n{H/a) 



1 



n' 



din a ^ -H2 

which, on using the background equations (3.3)-(3.4), and in order to get a exponential expansion phase, imposes 

«2 

|p + P|<l. 



(This generalizes the slow-roll condition for the scalar field in /c-essence model 



(3.14) 



(3.15) 



2-H^^^X 



<1.) 



4 Linear perturbations 



We now derive the second order action governing the dynamics of Hnear perturbations around the background solution for the 
general action given in (2. 1), using standard techniques (and in particular following [40]). One should note that this second order 
action has already been studied in [23]: taking into account differences of notation, we have checked that all our results agree 
with the corresponding expressions given in [23]. 

Before proceeding, it is useful to count the number of dynamical degree(s) of freedom in our system. Since the Lagrangian 
(2.1) is constructed such that the corresponding field equations of motion (for both the scalar field and the metric) are second- 
order, there are no additional degree of freedom: the dynamical degrees of freedom of the model are exactly those of standard 
General Relativity with minimally coupled single scalar field, namely one propagating scalar mode and two tensor modes. 

Given that the curvature perturbation C is conserved on large scales [39], we have found it more convenient to carry out 
perturbation calculations in the uniform scalar field gauge with 5(t> = 0. Then the perturbed metric takes the form' 



ds^ = a2(7y) [- (e^" - e-^'^S'^dil3djl3) dr]^ + 2d^l3dridx' + e^'^S.jdx'dx^] , 
so that in terms of ADM variables, the lapse, shift, and metric on spacial slices are given by 

N = ae°' ,Ni= a^d./S and hij == a^e^^(5„ 



(4.1) 



(4.2) 



respectively. Just as in minimally coupled fc-essence, the lapse and shift and thus a and (3 are constraints in our model (2.1). 
Thus in general a, (3 can be solved for perturbatively in terms of curvature perturbation (^: 



E 



a 



("')' 



/3 — 2^ P(m)i 



(4.3) 



where the subscript "(m)" denotes the order in C,. Furthermore, since a, /3 are constraints, their first order solution a(i) and /3(i) 
are adequate for our purpose to evaluate the power spectrum and bispectrum. 

4.1 Power spectra 

4.1.1 Scalar perturbations 

On using the form of the perturbed metric (4. 1), a straightforward Taylor expansion of the full action (2.1) together with some 
integration-by-parts and use of the background equations, yields the quadratic action for Q, a and /?: 



'S'(2) [C, ", /3] 



drjd xa 



-35cC + cc {dO - Sn^rua^a^ + 2g^da ■ dC + GHfaaC + ^gc^C'd^p - 2Hfo.ad^P 



(4.4) 



In (4.4) four dimensionless coefficients appear (that is, dimensionless when we reinstate the Planck mass A/pi), and they are given 
by: 



12 






1 - 4XG(') + 2G(2) - '4x ia'G'^':} - HG 



1 + 2G(2) + ^X 



:(3)w^' 



G 



['J ir ■ 



2x^(3) 



a 



6W 



K + AX^K 



XX 



W) + a^G\l 



3-^2 V ^^'^ 



G 



(i; 



n 



X[iG 



(1) 

X 



2XG 



(1) 



XX 



-G*^' - 2X (sG^y' + 2^ (7G|xx 



2XG 



(2) 

,xxx 



n 



4^2^(2) 
4-^ ^,XX4 



8XG 



(2) 
,X<p ' 



G 



(2) 



-?^x(25g5 



AX 6G 



.(3) 

,xx 



XG 



(3) 
XXX 



-6X{AX^G%, 



16XG 



(3) 

X<t 



9G 



(3) 



fa 



1 + |xgW + 2 (g(2) - AX (g|1) + XG% 



|(2XGg, + G|? 



GHcf)' 



.(3) 



.(3) 
,XX 



12X{2XG 



(3) 
Xcf, ' 



3G 



(3) 



Here 5^ etc are normalized such that in fc-essence models (K 7^ 0, G^^'^'^^ = 0), 



(4.5) 
(4.6) 
(4.7) 
(4.8) 
(4.9) 
(4.10) 
(4.11) 



5C = cc = /a 



TOq, = 1 — 



3-^2 ^ -^ 



and 

+ 2XK 



^^^"^ W 



'Throughout this paper, spatial indices are raised and lowered and summarized by 5ij and 9^/3 = &^^ didjP, 5^^ , and (9i/3)^ = S'^difidjfi, (didj/S)^ 



where c^ is the standard expression in fc-essence models, namely c^ = K -^j {K x + '^^^ xx) which we will derive below. 
Notice more generally that neither K nor G'^^^ contribute to g^ and cq\ thus in the G-inflation limit we also have gc, = c^; = 1. On 
the other hand G^^-' and G^'^-' contribute to all four expressions. Finally, in the standard Galileon limit ( 1 .4) the above expressions 
simplify a great deal. 

The constraints a and /? are now obtained by varying (4.4), yielding 






C, (4.12) 



.9C f,2. o / 5c™- 



a^/3=:^9^C-3(^^-ljC'. (4.13) 

(Again one can verify that in fc-essence limit, (4.12)-(4.13) reduces to the standard results [41] a = ^ and d^(3 ~ ^^f + 
^X {K X + '^XK xx) C = -^ + p-C'-) On substituting the constraints (4.12)-(4.13) solution into (4.4), we obtain the 
final quadratic action for C,, 

5(2) [C] = J dvd'xa''^ (C' - clidCr) , (4.14) 

which is controlled by two key coefficients 

which determine the normalization and propagation speed of ( respectively. In the fc-inflation limit, es = e whereas c^ is given 
by the familiar expression [9]. In G-inflation, our expressions agree with those of [20, 33] (on using the background equations of 
motion), while in the pure Galileon limit they agree with those given [28]. Finally, similar expressions for the general Galileon 
model can be found in [19]. 

Finally, in order to avoid ghost and instabilities, we must impose 

e, >0, cl>0. (4.17) 

Notice, however, that since e^ is not simply related to the slow-roll parameter e, requiring £„ > does not necessarily impose 
H <0. 

4.1.2 Tensor perturbations 

For completeness, we also consider the tensor perturbation 

ds^ = a^iri) [-dif + (% + h,j) dx'dx^] , (4.18) 

where the symmetric tensor hij is transverse and traceless. The quadratic action for the tensor modes hij can be easily derived 
as 



/'2 



(4.19) 



where 

eT = c^, 4 = ^. (4.20) 

5C 

Notice that no new parameters enter into the tensor perturbations, and that in both fc-inflation and G-inflation where c,^ = g^ = 1, 

we have er = ct = 1. 

In the general model (2.1), both the amplitude and the propagating speed of gravitational waves are modified relative to 

fc-essence model. In particular, the tensor perturbation can be either amplified or suppressed by tuning the parameters 5^ and c^^, 

which can be done only with G^^'^^ . Finally, in order to avoid ghost and instabilities, we require 

ffC > 0, cc > 0. (4.21) 

We now introduce the polarization decomposition 

h,jiv,k)^ Y, h-\,j,k)et^{k), (4.22) 



where the symmetric tensors e^'^ are transverse and traceless k^elAk) — ef^ — 0, and satisfy the orthogonal and normaHzation 
condition: efjk)efj*{k) = 26ss'- In terms of h", action (4.19) can be rewritten as 

5(2)T= E fdvd'x^'-^ihf^cUdhs)'). (4.23) 

4.1.3 Quantization 

We now quantize the perturbations following the standard procedure e.g. [42]. 

To this end, we assume a quasi-de Sitter inflationary background, and that the Hubble parameter H, eg, Cs, er and ct 
are slowly varying with time during the period when cosmological perturbations, whose scales are of the current observational 
interest, are generated and exit the Hubble scales. In this case, we use the familiar mathematical trick, i.e. first evaluating the 
power spectra by treating H, e^, Cs, et and ct are exact constant, then taking into account their time-dependence by identifying 
H = H{rj^,{k)), Eg = es{ri^,{k)) etc with —Cgkrj^, = 1 in order to evaluating the spectrum indices, up to the first order in 
slow-varying parameters. 

The canonically normalized variable corresponding to ^ is defined by (^ = —^/2e^( = z(, and from (4.14) it satisfies the 
familiar equation of motion 

Ck + (elk' --)a = 0. (4.24) 



On selecting positive frequency solutions which correspond to the standard Bunch-Davis vacuum deep inside the Hubble scale, 
and on converting back to (, we find 

a(ry) ^ — ^^ (1 + ^cskv) e^ ''=^ (4.25) 

(where H = H/a) so that the power spectrum on super-Hubble scales (Icsfcyyl ^ 1) is given by 

n = ^2 ' (4.26) 

where all quantities are evaluated around the Hubble exit at Icsfc??! = 1. Notice that this differs from the usual single field result 
in that e is replaced by e^ which, contrary to fc-inflation, may not necessarily be small (a larger e^ would tend to decrease the 
power spectrum). Similarly, the tensor power spectrum is given by 

2 zj2 

Vt = — , (4.27) 

which can again be amplified, or reduced, depending on the values of etCt- At this point, we emphasize that our expressions 
(4.26) and (4.27) are consistent with the corresponding results got in [23]. 

The tensor-to-scalar ratio therefore reads 

Vt 16esCs 



V(; exCT 

The spectral indices for scalar and tensor perturbations are given by 



(4.28) 



where in the above. 



71^-1= -^ = -2e - f^s - s, (4.29) 

d In fe 

dlnVr 
fiT = ,. . = -2e - j-jT - ST, (4.30) 

dliifc 

d hi Cs d hi ey d In Cs d In ct 

Vs = -7] , ^T ^ -r, , « = Tri ' ^T = -T--, • (4-31) 

a In a d In a d In a d m a 

2 



In terms of these slow-varying parameters, eg defined in (4. 15) can be recast as 

1 e^ 

£., = -p ^ (e + 1 + a^ (2?7T - sj - 4st)) - er, (4.32) 

where 

Sf = din fa/dlna. (4.33) 



Note that we have not managed to ehminate fa and use it here as a slowly varying parameter. 



In fc-inflation and G-inflaton, et = ct — 1, tjt — st — 0, ut — — 2e, so that the tensor-to-scalar ratio r = WegCg. In fc-inflation 
where Cs = e, then r = —ScsTit which is the standard consistency relation. In G-inflation e,, and e are related in a complicated 
manner, but one can again derive a modified consistency relation. This was already found in [18, 19, 23, 26, 30], and in our case 
is expressed as 



r = -8c. 



JaCrp 



riT 






- 1 



rjT + St ~ 2a^ (2'7t ^ Sf — 4st) 



(4.34) 



Thus, this gives the first distinctive feature of generalized Galileon model different from fc-essence model. 

5 Bispectrum of curvature perturbation ( 

5.1 Cubic action for the curvature perturbation 

In this section we again follow the standard approach and extend the previous calculation to third order in perturbation theory. 

Due to the presence of various higher-order derivative terms in generalized Galileon model (2.1), the number of cubic order 
terms in perturbative Lagrangian is very large! In order to group various terms, we have proceeded by the following steps. First, 
we integrate by parts to remove the time derivatives of a and /3, and second we simply (or are able even to eliminate) terms using 
the background equations of motion (see Appendix C). After rather tedious and cumbersome multiple integrations-by-parts, we 
eventually arrive at the following third-order action for C,, a and (3: 



•So) [C, ", /?] 



drjd X a I qq 



-9CC" + 2C' (Ca'/3 + 9,C5'/3) - a {d,Cf + {d^fif d\ - ^C {"^ad^ - {d'/sY + {d.d./Sy 



+ccC (9,0" - m^m^a-C + 2H/„« (9CC' - Cd^P - d^Cd'/3) 

n 

+X2a 



^'3 „ ^'2g2p ^ 1^, 1^4^52^ ^ (52^) 2 _ (^g^Q^f^f^^ _ „ ^g2^g2f^ _ 9,3, C^,;^, /3) 



This^ contains 28 cubic interaction terms and nine dimensionless coefficients: four of these are g^, c^, /„ and nia which already 



^We have compai'ed our result with (e.g.) the coiTesponding results recently presented in [26]. It is interesting to find that, (5.1) has essentially the same 
structure with eq.(46) in [26], which also has 28 cubic terms controlled by nine independent coefficients. More interestingly, these coefficients have exactly the 
same relations among them. 
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appeared in the second order action (4.4); and the remainder we call A^, {i — 1,- ■ ■ , 5), whose explicit expressions are: 
Ai 






A, = 1 



G^^^-4Xid^^+XG% 



i^X {5G^^ + 2XG%) - 12X {2XG% + 3G|^' 



A3 = 1 + j^X (sgJ + 2XG[]^x 



^(2) 



2 


G(2) - 


6H0' r 
a2 ^ [ 



.(2) 



^(2) 



2X ( 5G7 + 2X 7G)^x + 2^G ^^^^ 






4-'^^G,xx</. + 8-'^'^,x</ 



'^S') 



2bG'V + 4X ( 6G*iV + ^G*^^ 



:(3) 



.XX 



XXX 



12X 4X XG 



A4 = 1 + 2 (g(2) + 4X2g|^'^) - i^^ (3G^;^ + 2XG 



^(3) 



(3) 
XX 



,XX0 " 



12X 2XG 



4G 



(3) 



(3) 
X0 



9G 



(3) 



G?) 



A., = 1 



20' 



X 



[K - 2X {K^x + &XK^x^ + 4^2;^^^^;,)] 

2a2 



^(1) 



.(1) 

.XX 



.(1) 

.XXX 



^X[AX^G^^^^ 
3-^2 V .XX, 



8XG 



(1) 

,X4 



G(2) - 4X (TG^i' + 4X ( lOG^i'^ + 7X&^^^^ + X^G%^^ 






,XXX4 



G 



(2) 



+ £^xfl25G(!.)+218XG(^) 



-12X ( sx^G^^^^ + eox^Gg^, 



^2^(3) 

.XXX 



-3^(3) 

,xxxx 



98XG 



(3) 

,X4 



27G 



(3) 



G 



(1) 



(5.2) 
(5.3) 



(5.4) 
(5.5) 



(5.6) 



where again, Ai are properly normalized. 

At this point, it is useful to recall the fc-essence model values of A^'s, which are given by 



Ai = 0, 
As - 1 



sn- 



A2 — A3 — A4 — 1, 

[K - 2X {K,x + &XK^^x + ^X^K^^xx)] 



(5.7) 
(5.8) 



Using the background equations of motion, in fc-inflation models, A5 can alternatively be written as A5 = — ^^y A, where 

2 



\ = X^K 



XX 



-X^K 



XXX' 



which is the "popular" combination which was introduced (e.g.) in [41] (note the original definition of A has dimension as H^). 
Thus A5 is the natural generalization of A in the Galileon model (2.1), whereas Ai, A2, A3 and A4 are new parameters in Galileon 
model, which are trivial in /c-essence models. It is interesting to note that in G-inflation, Ai = 0, A2 = A4 = 1 as in fc-essence 
model. Their corrections arise only when higher order Galileon terms, i.e. G^^^ etc are included. 

Finally, we can eliminate a and /3 in (5.1) using the constraint solutions (4.12)-(4. 13). After another set of cumbersome 
integration-by-parts, the number of cubic interaction terms surprisingly from 28 to 10 and we obtain 

Ai 



n 



Ay 






(5.9) 



with 



- a-2w 



V' = d-X 



(5.10) 
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The 8 dimensionless parameters in (5.9) are given by 

„2X ;^ 



Ai 



ere. 1 , _^ / ^ ^sCr 



CtCJ fa ' "^ V" <^TC1 



/< 



"3'-l-2'4 



A3 et / e. 



f 2 „2 \ : 



,eT 



1 ej. A5 



A2 = 4 



1 ex 



(e + 77t - s/ - 2st - 3 - 77s) 



A3 

A4 






er 



— (1 + 77T + ?7s + e ~ s/ - 2st - 2s) + Cj./q ^ C^Vt + e- Sf + 1- Ast) 



3Ai 



^T A3 CT ( r,A2 A4 

7^72 + 72" 1 ^7 f 



'" 2eTrf 



A, 



1 + Ai^ (3 - si - 4sT - e + 277t) + -^ (3 - 2st - S2 - e + ?;t + s/) 



Afi 



A7 



4eTc| 



64/^ 



3 + Ai ^ (e + si + 4sT - 27?T - 3) + -^ (e + 2st + 53 - 77T - s/ - 3) 



1 + 3Ai-^ 277T + Si + 3e - 1 - 2s/ - 4st - 2 

A2 / ^s^l^J a 

+-7- I 3?7r + S2 + 3e - 1 - 6st - 3s/ - 3 g — 

/a V ^T 



As = 






er 



1 + Ai-^ si + 2e-2-s/-/„ 



-S^X 






1X0/ EcC 

"o 7^ 'yr + S2 + 2e - 2 - 2st - 2s/ - 2/„^ 



(5.11) 
(5.12) 
(5.13) 
(5.14) 
(5.15) 
(5.16) 

(5.17) 

(5.18) 



where e^, c^, ey and ct are defined in (4.15), (4.16) and (4.20) respectively, and the slowly varying parameters ris,r]T, s and st 
in (4.31) and (4.33). Finally 

.. = ^. (5.19) 

dm a 

Equations (5. ll)-(5. 18) exactly re-produce the corresponding values of A^'s in the fc-essence Umit, which are (recall in fc-essence 

limit, Cs ^ e, ex = Ct = fa = A2 = A3 = A4 = 1, Ai = 0, A,5 is given in (5.8) and is usually written as A5 = ~jp) 



A _ « n i\ 2A 



A3 = — (1 + e + 77e - 2s - c^) , 

A4 = A7 = Ag = 0. 



A2 = ^ 



3 + — (e - 3 - 77e) 



A5 = -:t:4(4-^): 



^'= = 4d' 



where ij^ 



din 



din 



In [41], the contribution to the bispectrum of curvature perturbation up to the second-order in slow-varying 



parameters were calculated, which correspond to our Ai, A2, A3 and A5, whereas Ar was neglected in [41] as it is third order in 
e. 

In the last term in (5.9), we introduced 



SC2 






which is proportional to the linear equation of motion for C obtained by varying (4. 14) with respect to (, and 

FiO = /iCC + ./2 [id.,cf - d-'d,d, {dXd^C)] + /3 [a.cav - d-^d,d, (9'ca»] , 

with 



/i = "^ 

^' = 4^5^Z^(2Ai/„ + A2|J, , 

/3 - ~o_2o,^ -272 Al/a + -2"^:: 



2amfa e^j,c 



(5.20) 

(5.21) 

(5.22) 
(5.23) 
(5.24) 
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5.2 Bispectrum of the curvature perturbation 

Having derived the full third order perturbation action (5.9), it is a standard but subtle'* exercise [4] to evaluate the corresponding 
bispectrum, i.e. the three point function of curvature perturbation C. The tree-level contributions to the three-point functions of ( 
from cubic interactions described by (5.9) can be evaluated using 



{CkAv)CkAv)Cksm = -2^ 



I I 'd7j'{o\aAv)CkAv)CkAv)Hdv')\o) 



(5.25) 



where in the right-hand-side of (5.25) all quantities are in the interaction-picture, jy; is some initial time when perturbation modes 
are deep inside the Hubble radius, which can be conveniently chosen to be rji = — oo. At cubic level, the Hamiltonian density is 
simply given by H(3) = -^(3)- 

In the following, we assume that the eight dimensionless coefficients Aj (i = 1 , • • • , 8) are approximately constant when 
evaluating the three point functions using (5.25). We simply collect the final results, and on denoting 



{Ck^CkAks) ^{2Trf6Aki+k2 + k3)BAki,k2,h) 
the bispectrum B^{ki, k2, k^) is given by 



(5.26) 



Br 



htO htO htO ^ 
l\i-\ fXicy A,Q Cg 



(fclfc2fc3)^ 



ai+|a4)+§ 25:^,^^1-1 5:fc.^fc3 



i>] 



i^'j 



+^ f E ^? + ^ E fc^^l - Ji E ^?^.^) + T f E fc? 4 E ^^^1 - Ji E ^f^ 



Ae 
'2if2 



2 E '''^ + E ^*^J " 3 E ^»^^J - 2fcifc2A:3 J2 ^^^1 



j#j 



i#J 



" \ J J>i / \ i>3 



3fciM:3 



^Ag 



-f^^ I '^^Y.kf + 3fcifc2fc3 E ^»' - 2 ^ fcf - 5kik2k3K^ - 12 51 ^^^1 



(5.27) 



where K = ki + k2 + k^. 

From (5.27), various contributions to the bispectrum Bq from different three-point interaction terms in (5.9) can be deter- 
mined. Notice that ('^ and C'^^^C contribute the same momentum shape. This can be understood as follows: since we are using 
the linear solution of ( to evaluate the three-point function, we can replace d^C in C'^d'^C, using the linear equation (5.20), i.e. 
c^gd^C, = 2T-LC,' + C" (note we are approximating H, Cs, Cs, A^'s etc as constant). Thus we can write 



'^' • a/2q2/- _ '^' a /-I2 1 



a2 2A 



4 A/3 



_A.c"#C = 5pA,c«jj(2«C' + C») = -^r, 



which gives the correct coefficient =%*■ in the first term in (5.27). In fact, for a similar reason, the shapes proportional to A7 § 

appearing in (5.27) are not new and must be a linear combination of the others^. Indeed, on using the linear equation of motion 
one can also integrate by parts the last two terms in (5.9) which become 



As 



d\id,0 -CdrdAd.CdjC) 



a'^ [d\diCdi4> - Cd,d, (d.Cdji;)] 






n ci "^ ^^ 



1 



1 



1 



-Ag-CC" - xAsC (5.0 + ^SAsC'S.CS.V 



c^ 



'^We have only calculated the bulk contributions to the third order action (5.9) and not kept track of the boundary terms. This should give the correct result 
for the bispectrum: see [4] and [5] for a discussion of this point in the case of fc-inflation. 
^See [44]. We thank S.Renaux-Petel pointing out this important point to us. 
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/A4 


+ 


DVV7 


9A7 


+ 


BAg 
2cr 


3A7 




As 


c2 




2 ' 


3A8 
^2 


■} 


A, 



Effectively this means that we can set A4, A7 and Ag to zero in (5.27) and replace Ai g by 

Ai = Ai 
A2 = A2 
A3 = A3 

A5 = A5 + ^^ , Afi = Afi 
Finally, one may be interested in evaluating the dimensionless non-linear parameters, defined by [43] 

/NL(fcl,fc2,fc3) = 1 ^^T^ ' . (5.28) 

6(2.2)2^2^3^^^3 

To get some idea about the amplitude of the bispectrum, here we choose the equilateral momentum configuration fci = k2 = k^ 
and obtain 

5c?/, ,, 51 A3 2, ^, „, 39 A7 SlAgN ,, ^^, 

= ^^(A,+6A2 + -^ + ^A4-3A,-3Ae--^--^j, (5-29) 

where A^'s are given in (5.11)-(5.18). 

6 Conclusion 

In this paper we have studied cosmological perturbation theory in the "generalized Galileon" model proposed in [6, 7]. In 4 
space-time dimensions this is the most general scalar field theory whose Lagrangian contains derivatives up to second order 
while keeping the equations of motion which are second order and lower. Our model, which includes fc-inflation as a special 
case, represents a large class of scalar field models which has not been investigated so far 

The present work is the first step in exploring the cosmological implications of "generalized Galileon" models. By deter- 
mining the most generic second and third order actions for the curvature perturbation, we calculated the power spectra of scalar 
and tensor perturbations as well as the bispectrum of the curvature perturbation. We found modifications of both the ampli- 
tude and propagation speed of tensor perturbations due to the presence of £'^^' and £('^^ Galileon terms. Correspondingly, the 
tensor-to-scalar ratio is modified relative to fc-essence (see also [23]). We have also showed that, although there are higher-order 
derivatives in the Galileon model relative to /c-inflation, there are no new contributions to the bispectrum of the curvature pertur- 
bation. In order to get a feel for the strength of the bispectrum, we evaluated the non-linear parameter Z^^" for the equilateral 
configuration ki = k2 ~ k^. 

In the future, the strength and shape of primordial non-gaussianities will be constrained both by CMB data as well as, 
independently, by galaxy clustering data. Indeed, it is expected [32] that upcoming observations of high redshift clusters may 
possibly enable one to put limits on /nl at the level of a few tens. It remains a very interesting and open question to see whether, 
as a result of this different data, strong constraints may be put on the different unknown functions appearing in the generalized 
Galileon model. Our result (5.27) will be a starting point for such an analysis. 
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A Equations of motion 

The explicitly second order equations of motion for (p are obtained from (2. 13), and are given by 



n=0 



(") 



0, 



with 



£(0) 


= acf>K\x + K^xx^^.XV^4> - 


-2XK^x^ + K^^, 


£{l) 


= 2n0Gy - 2xg;^^ + g5 


"(00)2 - (V^V,(/.) 2 - i^^.V^f^V^f/)" 




-2G^^^, {XU^ - V^XV^0) + G^i'^ (V^XV^X + U^V,,XV^i 



(A.l) 



(A.2) 



(A.3) 



and 



£•(2) ^ G 



(2) 



(D^)^ - 300 (V^V,0)2 + 2 (V^V 



/J- V 2> 



2i?p^a.V''0V"0V^V'' 



-G 



(2) 



2n0v^xv^x + ( (00)2 - (v,,v,0)2) v,xv^ 



-2G'^^G^,V^V^ 



2G 



RG 



(2) 



G' 



(2) 



(<^ 



3 (Dg 



(2) 

,XX0 

2 



2V^XV'-'XV.V^<i 



^ (V^V^0) ^ - (00)^ + 2V^XV^X + 2n0VpXV^(; 



(V^V^<; 



2i?X - 4i?,,^V^0V% 



^(2 



2G^x'^. (VpXV^0 - 2X00) . (A.4) 



Though it is not particularly eluminating, the expression for S'-'^^ is given by 



£•(3) ^FiG''^^+F2G 



.(3) 

,xx 



FsG 



(3) 



3"-^. XXX 



FaG 



(3) 



4'-^,XX<i 






i^sG 



(3) 



where Fi , • • • , Fg are scalar functions given by 
Fi = 3R 



(00)2 - (VpV,0)2 - ii'^,V^0V^ 



-12i?,,, (V''V'^0VpV 



3 (2i?ppi?^ + 2RP^R^p,, - R/''^R,pax) V'^0V% 



O0V'^V''0) + Gi^^p.^V^V^VV" 



F2 = 3i?(VpXV^X + O0VpXV^ 



i?u.V^0V^ 



Fs 



- 3 ^(00)" - (VpV,, 

-6i?p^ (2V''0V''Vp0V''X - V^XVX - V^V''0VpXV''0 - 2O0V^XV''0) 
-6i?pp^^ [(O0V^V''0 - V^V^0V''Va0) ^0^0 + 2VP(j)V''XVV4>] , 



(00)2 - (VpV,0)2 



(O0)^ - 300 (V^V,0)2 + 2 (V^V,0)' 



F4 = 12 (O0.gp, 

-2X ' 
F5 =4 



V„Vi.0) V^XVX + 6 



(O0)'-3O0(V^V,0)2 



VpXV.; 



(O0)2 - (V^V,0)2 



VpXVc? 



p V ,^, -r 2 (VpV,0)' 

6i? (XO0 - V^XVp0) 
l2Rp^ [XV^V''0 - O0V^0V''0 - 2V^0V''X] + 12i?p^apV''0V"0V^V'^(; 



(O0)-'' - 3O0 (VpV,0)2 + 2 (VpV,0)^' 



Fe := -6X 

F7 : 
Fs : 



(O0)" - (VpV,<; 

12Gp,V^V>, 
6Gp,V^0V'^0. 



eVpXV'x + 6O0VpXv^ 



(A.5) 



(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A. 10) 
(A.ll) 
(A.12) 
(A.13) 
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B Explicit expression for T^u) 

The energy-momentum tensor corresponding to C'^^^ takes the following form: 



l^iU 1 



where Ci , • • • , Cio are scalar coefficients which are given by 



Ci = 6G 



(3) 



(VpV 



p V o-V 



-6G 



-3G 



(3) 

(3) 
X 



{U(j}Y + 26^,^0^^ -RX - GC^Z (VpXV^ - 2xn,^) 



(3) 



X (VpV^0) ^ - {n(j)y + 2VpXVPx + 2n(/)VpXv^ 



(n,j 



-47?^^^^^ 



3n(/)(VpV„ 



i^VpXV^ 



2(VpV,^ 



2DcbRp„yPcj)V'' 



2i?pAarV''(/.V"0V^V^ 



-G 



(3) 



3 ( (00)^ - (VpV„<?i) 2 ) VaXv^0 - evxvxVpV^;/) + en^v^xVpX 



(B.l) 



(B.2) 



.(3) 



Co = 6G 



(3) 



{Ocbf - (VpV 



p V (jV 



.(3), 



:(3)i 



6G^;^i? + 6G^;;n0 - gg^^^Gp^vp^^ 



+G 



(3) 



(00)3 - 3 (VpV,0)' 00 + 2 (VpV„ 



C3 = 12G^3^ 



6G^i^^O0, 



GG^'^^i?- 



6G 



.xx"- 

(3) 



XX 



(O0)2 - (VpV,^ 



24Gg^O0- 



12G 



(3) 



.(3) 



G5 = 12GT0(i- 12G^7.^ - 6G 



(3) 



.(3) 



00" 



.X 



eVpXVc; 



Ge 

G7 
Gs 

Gg 

Gio 



^2g5^ (6X0,? 
-12G(^)^O0-24G(^)^, 



(VpV,0)' - (00)' + Rp^VPcj^V'i 
f 6G^^'x (VpXVX + DcpVpXyPc 



12G 



(3) 
XX' 



-UG'-j^Dq 



^2G%^- 



-(3) 



(3) 



12G^^^ - 6G^^^xV„XV>, 



42G 



(3) 



XX' 



12G 



(3) 
X' 



and r/,,y represents terms whose fi, i^-indices explicitly depend on or coupled to Riemannian tensors: 



r(3) 



-6G 



.(3) 

,0 

(3) 



12G^'' {XRp, + Rpp^^y^V' 



' 2V(^0i?,)pV''0) 
V^V^vVi^Ap.CM - i?p. V''0V'^V(p 



{2[(O0V''0 + V''X)i?p(p 
~VpXVPcj)R^, + 2V''0i?p(pV,)X 
-i?p(^.). ^0 (O0V^0 + 2V^X) - 2V>VVi?A.p(MV"V,)0}. 



V. 



(B.3) 
(B.4) 
(B.5) 

(B.6) 
(B.7) 
(B.8) 
(B.9) 
(B.IO) 
(B.ll) 



(B.12) 



C General consideration of the structure of perturbative action 

Due to the complexity of the perturbative calculation, it is useful try to deduce some general properties of the perturbed action. 
As we now show, on using the background equations of motion, two types of cancellation occur rather generally. As a result 
the conclusions of this appendix are useful not only to simplify expressions but also provide us consistency check in practical 
calculations. 

We first argue that in the quadratic Lagrangian there are no "mass term" C,"^ and aC,', and similarly in the cubic Lagrangian 
there are no <^^ and a^^. To this end, we focus on the corresponding fully non-perturbative Lagrangian of (2. 1), substitute (4.1), 
and furthermore neglect all spatial gradient terms. This yields 



C « e3« [A{C, a) + a'BiC, a) + C"G(C', a)] 



(C.l) 
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with 



A= —e 



a e 



B 



C 



a^e^'^K + QG\'^{n + C') (00 ((2H - C,'W - Sf/-") + 

■*" (sCH + (,'){% + 2C') + iW + 6G(2) (ch + C){n + 2c') + n') - g^^^ {{2% + k'W + 0")) 

6a2e2"(H + C') (-G JV (C'0' + 0") + SG^^) ((3C'(H + C') + 2H') 0' + (H + C')0")^ 
(H + C')' (0')' - a^e^" (3 (1 + 2G(2)) (H + C') - gI^V) 
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(^)( 



+6a2e2"(H + C')0' (9G(3)(H + C') - G^^V' 
e-3" (3a2g2a 1^^ ^ 2G(2)) - 36G(^)(^ + C')0' 



(C.2) 

(C.3) 
(C.4) 



We emphasize that in (C.2)-(C.4) we neglect all spatial derivatives in the functions K and G'^^'^''^\ so for example K 
K{^ '^^2 4''^, 0)- Notice that ( enters the Lagrangian (C.l) only through an overall prefactor e'^^. 



Although (C.l) can further simplified through integration-by-parts, it is already adequate for the following discussion. The 
perturbative expansion of (C.l) yields 



£ = £ 



(0) 



^(1) + ^(2) +-C 



(3) 



where the linear action is 



£(1) = (3^ - A[^, + C") C + (A,„ - B') a, 
which is vanishing and yields the corresponding background equations of motion 

£^ = 3A- A[^, + C" = 0, 



(C.5) 
(C.6) 

(C.7) 
(C.8) 



where quantities are understood as their background values. After integration-by-parts, in the quadratic Lagrangian terms which 
which are proportional to C^ and (a are 



and in the cubic Lagrangian: 



£(2) ^::{5A- A[^, + G") C' + 3 (A,„ - B') (a, 



£(3) 3 I (3A - A[^, + G") C^ + I (A,„ - B') ('a. 



(C.9) 



(CIO) 



2 ^ ■'> ' ' 2 

At this point, it is explicit that C^ and a(^ in the quadratic Lagrangian (C.9), C^ and aC^ in the cubic Lagrangian (C.IO) are 
proportional to the background equations (C.7)-(C.8) respectively, which are thus vanishing. 

The above analysis can be directly generalized to the higher orders. Generally, in n-th order perturbative action, C" the 
proportional to the background equation of ( (C.7), aC"~^ is proportional to the background equation of a (C.8), which is the 
background energy constraint. Thus both C," and aC"^^ must be vanishing in the n-th order perturbative action^. 

Another type of cancelations is the self-interaction terms of f3. One can explicitly check that the coefficient before /3^ 
self-interaction terms is proportional to the background energy constraint and thus vanishes. 

We wish to emphasize that the above cancellation is a virtue of choosing ADM-compatible perturbation variables^, i.e. 
N ~ ae" and N = a?diP. If one alternatively chooses other variables, in general there is no such cancellation. 
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